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Abstract 


We derive the odd parity perturbation equation in scalar-tensor theories with a non minimal 
kinetic coupling sector of the general Horndeski theory, where the kinetic term is coupled to the 
metric and the Einstein tensor. We derive the potential of the perturbation, by identifying a master 
function and switching to tortoise coordinates. We then prove the mode stability under linear odd- 
parity perturbations of hairy black holes in this sector of Horndeski theory, when a cosmological 
constant term in the action is included. Finally, we comment on the existence of slowly rotating 
black hole solutions in this setup and discuss their implications on the physics of compact objects 
configurations, such as neutron stars, 
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I. INTRODUCTION 


Black holes are among the most interesting and fascinating objects in gravitational 
physics. They form by gravitational collapse of heavy dying stars, but the thermodynamical 
properties they display makes contact to some fundamental aspect, beyond the classic held 
theory. Indeed, black hole thermodynamics lies between the classical and quantum pictures 


of nature 


a. 


However even if black hole thermodynamics represents one of the most 


important achievements in order to describe gravitational phenomena from a quantum held 


theory point of view 


Q. B. 


a complete theory of quantum gravity is still unknown despite 


the numerous attempts in the last decades m (and references therein). 

Uniqueness theorems are well established for four dimensional asymptotically hat sta¬ 
tionary black holes in electro-vacuum. Israel’s theorems adressed the uniqueness of the 
Schwarzschild and Reissner-Nordstrom black holes ^1, Wald and Carter extended the result 

□ n 

to the rotating case [10, iUj. In other words four dimensional electro-vacuum stationary 
black holes in asymptotically hat spacetime (and regular in the domain of outer communi¬ 


cations) are characterised by only 3 parameters 


are described by the Kerr-Newman metric 


‘back holes are bald’ (they have no hair) 12| . 


12 


mass, angular momentum and charge), and 


l4. This result seeded the famous phrase 


15|- 


Appart from the asymptotically hat case, no dehnite result were found, except for few cases, 
for instance for the Kerr (A)dS family. Given a cosmological constant, this family of so¬ 
lutions depends on two parameters; the mass M and the angular momentum J, which are 
dehned asymptotically as conserved charges of the spacetime. These quantities obeys the ex¬ 
tremal condition \J\ < GM"^, whose saturation gives a degenerated horizon, while otherwise 
the spacetime is nakedly singular. 

Uniqueness theorems are not only important per se, but, have also inspired the famous 
Wheeler conjecture “Black holes have no hair” [^. This conjecture, established as a theorem 
in a large variety of cases, is one of the most important and studied theorems in black holes 
physics. Wheeler stated that gravitational collapse leads to a black hole conhguration in 
such a way that the hnal state can only be characterized by mass, charge and angular 
momentum. Any other additional parameters being regarded as “hair” are not allowed in 
the external hnal description of a black hole. 

The no hair conjecture was the arena of the beginning of a frenetic search for hairy black 
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hole solutions. Indeed, after the announcement of the conjecture the community started to 
look for solutions with different contents of matter with the aim of obtain some sort of hair. 


The hrst hairy solution 17|] was found in the context of Einstein-Yang-Mills theory. Many 
attempt using scalar helds in asymptotically flat spacetimes were done given origin to no 
lair theorems, as is the case for minimally coupled massive scalar helds and proca helds 

18|,|l9|. 


However, it was recently found in 


20 


2l| that hairy black holes in asympotically hat 


spacetime do exist, in the presence of a complex massive minimally coupled scalar held. In 


this case, the scalar hair is supported by the rotation, and no static limit exits. See 


22| for 


a review. Similar conclusions for complex spin-1 helds are suggested in [23|, due to the same 
mechanism than in the scalar case. 

Scalar tensor theories have been since then a fruitful theoretical laboratory for this re¬ 
spect. The hrst scalar hairy solution was found in 2^ and independantly in 25|] by consid¬ 
ering a conformally non minimally coupled scalar held with self interaction, unfortunately 
the scalar held conhguration is not regular at the horizon. Adding a cosmological constant 
this problem can be solved and other families of solutions have been studied 261-129 1. 

The most general scalar-tensor theory with second order equations of motion for the 
metric tensor and the scalar held has been introduced by Horndeski (dubbed Horndeski 
theory) 30j some decades ago. The model has been forgotten untill recently, where it has 
drawn a lot of attention, es pecially in the case of cosmology. In fact many inhationary 
models have been proposed |3ll437l| . The theory has been also investigated deeply in the 
context of hairy black hole solutions, with special emphasis on the case of non-minimal 
kinetic coupled scalar helds, particularly in the case where this coupling is given trough 
the Einstein tensor. Non-hair theorems for the existence of asymptotically hat black hole 
solutions have been presented [38|, l39|, nevertheless to date, many solutions to this system 


were found in the case of anti-de Sitter asymptotic behaviors 401-14Tj (and references therein). 
Compact conhguration have also been constructed in this sector, leading to constraints on 


the parameter space o 


It was shown in 


the model 481. 


49| | that the Horndeski Lagrangian can be written as 


L = p) - Gsicj), p)n0 + G4(0, p)R + G4,p(0, p) [(□(/>)' - (V^V.0)2] 

+ G5(0,p)G'm^^V^V> - ^[(D^)' - 3n(/.(V^V.(/.)2 + (V^V.0)3], (1) 
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where the functions K and Gi depend on the scalar held and its usual kinetic term p = 
This form of the Horndeski Lagrangian allows us to recognize more clearly impor¬ 
tant subsectors of the theory, like general relativity, Brans-Dicke, the non minimal kinetic 
coupling sector, etc. We will mainly focus on the subsector described by the following action 


S[g^lu, 0 ] = / 


R-‘2A 1 , ^ 


( 2 ) 


This paper is devoted to derive the perturbation equations for static spherically symmet¬ 
ric spacetimes, with the main purpose of determinate the mechanical stability of black holes 
in this sector under odd-parity perturbations. To this end we shall follow the very gen¬ 
eral technique displayed in 5^ by Chandrasekhar for the axial perturbations of spherically 
symmetric spacetimes. The stability of general relativistic black holes have been studied 
for a long time, starting with the pioneering works of Reggee and Wheeler Ml- Might 


they have been unstable, their astrophysical importance would probably have been depre¬ 
ciated. Instead they are stable towards linear perturbations. And even more importantly, 
they are characterised by quasinormal modes (see Ml for a review) because of the horizon’s 
absorbant nature. Quasinormal modes are eigen modes of the linear perturbation equa¬ 
tions with a complex frequency. The real part of the mode provides a timescale for the 
perturbation damping, while the imaginary part is the oscillation frequency. In particular, 
the quasinormal mode ringing describes very well the physics of a binary black hole ’after 
merger’. 

When considering black hole solutions in a realistic scenario and establishing their possible 
realizations in nature, the question of stability is essential. For an alternative model of 
gravity to be viable, conhgurations known to exist in nature should be well described. In 
this context, black holes are indeed thought to exist in nature, for instance supermassive 
black hole seem to play a very important role in galaxies formation, and solar mass black 
holes are suspected to be present in some binary systems. 

In the context of general scalar-tensor theories with second order equations of motion the 


group of Kobayashi et al. have made an extensive relativistic analysis of linear per 


Indeed they have considerer the full Lagrangian ([T]) and studied both, the odd [Ml] and even- 


urbations. 


parity 


54j | sectors of perturbations deriving perturbation equations from the second-order 


actions. Perturbating directly on the action allows them to formulate no-ghost conditions 
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for this kind of solutions. On the other hand in both works ^ and the authors also 
obtain general conditions, necessary but not sufficient, to ensure the gradient stability of 
static, spherically symmetric solutions. 


Note that unstable b 
vant in some context 


ack holes with a very long unstability timescale might still be rele- 


55[ |. otherwise, one could reasonably expect a viable alternative model 


of gravity to provide stable black holes. 

This paper is structured as follows. In section [TTl we detail the sector of Horndeski theory 
that we are interested in, and present the odd parity perturbation equations. In section 
mil we present the parity odd metric perturbation in this model. Next, in section IIVI we 
discuss the stability of some black holes solutions in this model. Finally, we discuss the 
particular case of the dipolar mode, i.e. the slowly rotating case in section |V] before giving 
our conclusions in section I VI I 

II. NON MINIMAL KINETIC COUPLING SECTOR 


We consider the following action, 

1 


S[gt^u, ^] = j 


{R - 2A) - - - V{cj>) 


(3) 


where is the Einstein tensor and a, /3 are positive constants to avoid the presence of 
ghost. The equations of motion for the metric and the scalar held (j) are 


and 


Rfiu G T A.g^i/ K T fiQ^u\ 0, 


{ag^'^VA - fdG^^'^VuA - = 0 , 

d(p 


(4) 


(5) 


respectively. In order to write down the equations of motion we have used that V^0, 

the compatibility of the metric Vp,g^^ = 0 and the conservation of the Einstein tensor 
provided by Bianchi identity V^G^'^ = 0. For simplicity in the notation of equation (|1]), we 
have dehned the following tensors 

( 6 ) 


2 a 




- (V^V»(V.V„0) + (V^v.0)n0 

+ gtiv 


+ ^(V“V^0)(V„V^0) 


(7) 
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III. PARITY ODD PERTURBATIONS 


In this section we analyze the odd-parity sector of pertnrbations aronnd spherically sym¬ 
metric and planar black hole solutions of the action ([3]) in asymptotically (anti-) de Sitter 
spacetimes. Once we have arrived to our master equation we will analyze the stability of the 
solutions making use of the Fourier decomposition of the master variable. Then our stabil¬ 
ity criterium will be based one the positivity of the spectrum of the Schrodinger operator, 
following the lines of 561. 

The perturbed metric reads 

dz^ 


ds^ = —A{r)dv + B{r)dr^ -|- C{r) 


-|- (1 — kz^){dip -f kidt + k2dr + k^dz)" 


, ( 8 ) 


1 — kz"^ 

where fci, ^2 and k^ are functions of {t,r,z). A{r), B{r) and C{r) are the metric functions 
parameterizing the most general static background solution of a scalar-tensor theory and the 
parameter k = ±1, 0 characterizes the topology of the spacetime. We consider the following 
ansatz for the scalar held 

4> = Mr) + (^^{t,r,z), (9) 

where 0o is the background held. Considering the Einstein held equations only at hrst order 
in e, we hnd that 

dt^ fd 


j^t Mo 
= en—— a- 
dr 

j^r Mo 

dr 


A 


AArCr ICl 

ABC\2 A ^4C Crdr]dt^ 


+ OM) = 0, 


( 10 ) 


a- 


B 


lArCr IM 1 
B^C \2 A ^ 4~C 2 


/3 


ACr + CAr 

A 


dr 


+ O(e^) = 0, 

( 11 ) 


these equations imply that 9*$ = 0 = dz4>. Substracting the equations E^ and EJ^, we obtain 




2a^a,<S> - P 

dr 


M Mo Cr ddcfpQ 2k d(j)o 


BC^ dr BC dr^ C dr 


3 BrCr d(j)Q Crr d(j)Q 3 CrAr d(j)Q Cr d(j)Q 


dr \ dr^ 


2 B^C dr BC dr 2 A^BC dr BC dr 
this equation imply that dr^ = 0, using these results in the equation E^ we hnd 

El = + 0{e^) = 0 , 


+ O{e^) = 0, (12) 


(13) 


where Vi arises from the expansion of the scalar held potential around the background 
conhguration, 

rlV 

(14) 




(t>=4>o 
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Thus, it follows that $ = 0, as obtained in the case of a minimally coupled real scalar field 


57l |. Using the zeroth-order equations, it is possible to check that the remaining equations 

are satished up to linear order in e if the following system is satished 

2 ' 


= — 
^ dz 


A 

C 


dt 


1 + 


1 + 

Pk f d(f)Q 


fdn ( d(j)o 


2B \ dr 
2 


il-kz) {d^k 2 - drks) 


2B V dr 


(1 - kz"^) {drki - dtk2) 


= 0 , 


(15) 


= — 
^ dz 




(3k ( d0o 


+ 


dr 




\/AB 


1 + 


1 - 
(3k [ d(j)Q 


2B V dr 
2 


{l-kz) {d^ki - dtks) 


2B V dr 


F^ = — 




(3k f d0o 


dt 


Cy 


1 - 


1 + 

(3k (d((>Q 


2B V dr 
2 


(1 - kz^){drki - dtk2) 


{d^k2 - drk^) 


= 0 , 


(16) 


{dtk^. - dM) 


= 0 . 


(17) 


2B \ dr 

As we can see, at hrst order in e, we do not have cosmological constant in the equations of 
motion. By introducing the variable 


Q = cJ 4^(+)(l - kz^f{d,k 2 - drks), 
V -D 


where 


the equations (IT^ - ffTT)) can be written as follows 

\/A dQ 

C^^B{l - kz^yP(_) dr 

The combination 9^(120]) + 9^([2T]) can be written in terms of Q, 


= -d^ks + dtd:,ki, 


^ d 


^^AB 


dr 


^ 1 dQ 


d 


+ {l-kzy^ 


dQ 


cy B P(_) dr 

Using the separation of variables Q = Q{r,t)D{z), we obtain 

1 [a 1 dQ 


dz [(1 — kz"^) dz \ A 


= ^d!Q. 


d 




Cy B P(_) dr 
1 dD^ 


-fQ = jd^Q, 


dz [(1 — kz'^) dz 


= —'jD. 


(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 

(24) 
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Let us focus on the case with k = 1. Setting z = cos 6* equation fl2^ allows to identify 
= D{z) with a Gegenbauer polynomial with 7 = (/ — 1)(Z + 2) where I > 1 must 
be satished. The Gegenbauer polynomials in terms of the Legendre polynomials read 


^- 3/2 

'-^i+2 


(^) 


sin^ 9 


d 1 dPi{9) 
d9 sin 9 d9 


_/2 

In this case we introduce the master variable = [GP(_)] Q{r,t) where 


4^. Using these dehnitions in equation fl2^ one gets the master equation 


_a 

dr 


( I d^C _ 

Qj.*2 

1 dC d7^(_) 
2CV(-) dr* dr* 


3 UC 
4G^ \dr* 


A 


+ 

T : 


1 d"P(-) 

2P(_) dr*"^ 


3 (dVi-)V 
4Pf_) V dr* ) 


(25) 


The mode stability can be explored using the Fourier decomposition of the master variable, 
^ = f which yields 


TfT., : = 




+ 7 


A 


1 d^V^.) 


CVm 

1 


+ 


dr* 


4G2 
dC dP() 


1 d^C 
2C dr*"^ 


+ 


471, 


'Tf’i-) 

dr* 


2P(_) dr 


*2 


2GP(_) dr* dr* 








where we have dehned It is important to note that even if the scalar held 

perturbation vanishes, equation fl26|) depends on the backreaction produced by the scalar 
held. If we choose /3 = 0, A = 1 — 2m/r and C = r^, this equation becomes the Regge- 
Wheeler equation. The spectrum of the operator "H is positively dehned as follows 


(26) 


dr*{^^)*n^^ = / dr* [I P +U 5 I P] - (^0.1^^.) {Boundary, (27) 


where 77 = ^ + 5 and 


dr* 


Vs = V + ^-S\ 


(28) 


If we choose S = we hnd 


Vs = 1 


A 


cn 


(+) 


(29) 
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IV. APPLICATION TO BLACK HOLE SOLUTIONS 


For spherically symmetric spacetimes the action ([3]) shows interesting black hole solutions 


27| 


k f a + (3A.\^ 

A{r) = — + -^/^[ -^ ^ 

a \a — pA) r 


^a.j3k 


r a- 13 k ^ ^ 


B{r) = 


a^((a — j3k) + 2f3ky 
[a — /3Ay{ar'^ + /3kyA{r) 


(31) 


(i 0 (r) 


{a + /3A) ((a — /3A) — 213k) 


(32) 


dr J k( 3 [a — I3A)^ [ar‘^ -\r I3k)^ A{r) 

where (3A ^ a, being /i a constant, LA := a/3/5 the (A)dS effective radius, C(r) = r and 
k = ±1 which characterizes the topology of the spacetime. The case /c = 0 integrate in a 
different manner given 

' d(j){r)p {a +/3A) 


dr 


Ka(3A{r) ’ 


where A{r) is simply 


* ' L 2 r B(r) 


(33) 


(34) 


In the expressions (13^ and (133|1 we must impose the condition, a + j3A < 0, in order to have 
a real scalar held conhguration outside the event horizon. 

For this solution, we obtain that V(+) > 0, hence because / > 1, 7 > 0 implies that V 5 > 0. 
Then in any region whenever A > 0, this implies that all the spherically symmetric four 
dimensional hairy conhgurations of this family of solutions are mode stable under odd-parity 
perturbations. In order to be true this statement we must require 


\Boundary 0, 


(35) 


i.e., vanishing perturbations at the horizon. 


V. SLOWLY ROTATING BLACK HOLES 


The aim of this section is to attack the existence of slowly rotating black hole solutions. 
In our analysis this can be archived by setting k 2 = k^ = 0 and ki = u{r) in equation (16). 
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Then the frame dragging function satishes 


dr 


L ^ / d(j)o ' 
VAB 1 ^ 2B\ dr 


.2^ d 




= 0 . 


(36) 


As concrete example let us take the case where k = 1 and z = cos(6'). Contrary to what 
happens for minimally coupled scalar helds ^ , the frame dragging equation also depends 


explicitly on the scalar held conhguration and not only through its backreaction. Inserting 
the solution (30)-(32) we surprisingly obtain for a;(r) 

C2 


a;(r) = ci + 

the same result than in standard Einstein gravity. Note that an interes 
analysis of this results has been recently addressed by Berti et al. in 


(37) 


ing and deeper 


58| . We comment 


however on the consequences of this result for the physics of compact objects, such as 
neutron stars. Since the frame dragging equation for compact objects is written around a 
static conhguration, in vacuum, the equation is strictly the same as the one we just presented 
above. As a consequence, the solution presented here, i.e. the same as in general relativity, 
is also the vacuum solution around a slowly rotating neutron star. As a consequence, this 
sector of Horndeski theory cannot be distinguished by observing (slowly rotating) neutron 


stars from a modihcation of the matter structure. This was discussed in details in 
non-rotating solutions, and this results strengthen the conclusions given there. 


48 1 for 


VI. CONCLUSIONS 


This work was focused on the analysis of odd-mode stability of black holes in the non- 
minimal derivative coupling sector of Horndeski theory. This sector is characterized by a 
non minimal coupling between the kinetic term of the scalar held degree of freedom and the 
Einstein tensor. For such a coupling and when the dynamics is governed by the linearized 
equations of motion, we showed that for a real scalar held conhguration and spherically 
symmetric spacetimes the solution described in is mode stable under odd-parity pertur¬ 
bations. It is important to note that, this can be stated because for our solution V(+) > 0 
always holds. This is due to the fact that our scalar conhguration is real provided by 
a -f /9A < 0. This rehects the importance of the inclusion of a negative cosmological con¬ 


stant. For instance in 


40| the solution does not contains any cosmological constant term 
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in the action, making the square of the scalar held derivative negative and in consequence 
P(+) > 0 is not always ensured. Such solution could then exhibits unstable behavior under 
odd-parity perturbations. In order to obtain our results we have followed the very general 


treatment given in 


50| . With a different approach, 


m 


53| . and without considering explic¬ 


itly the shape of the solution, conditions on the stability and no-ghost conditions in this 
kind of solutions have been presented before. 

Our result does not depends on the self-interaction term for the scalar held, behaving similar 
to what has been done for real minimally coupled scalar held in j^. 

The approach developed here allows the construction of slowly rotating solutions. Inter¬ 
esting enough we observe that the frame dragging function behaves exactly like in general 
relativity. This result have been explored deeply in the context of non hair theorems for 
galileon slowly rotating black holes in 


and has important consequences in the study of 


rotating compact conhgurations such neutron stars. Indeed, this type of solution can be used 
as the exterior solution of rotating neutron stars in the context of Horndeski theory and in 
this way search for new constraint on this model. The physical result can be anticipated here 
from the result we presented: the nonminimal kinetic coupling sector of Horndeski cannot 
be distinguished from general relativity outside a (slowly) rotating body. More details and 
applications to neutron stars will be presented elsewhere. 
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